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ABSTRACT 

Broad absorption lines (BALs) in quasar spectra identify high velocity outflows that 
likely exist in all quasars and could play a major role in feedback to galaxy evolution. 
The variability of BALs can help us understand the structure, evolution, and basic 
physical properties of the outflows. Here we report on our first results from an ongoing 
BAL monitoring campaign of a sample of 24 luminous quasars at redshifts 1.2<z<2.9, 
focusing on Civ A1549 BAL variability in two different time intervals: 4 to 9 months 
(short-term) and 3.8 to 7.7 years (long-term) in the quasar rest-frame. We find that 
39% (7/18) of the quasars varied in the short-term, whereas 65% (15/23) varied in 
the long-term, with a larger typical change in strength in the long-term data. The 
variability occurs typically in only portions of the BAL troughs. The components at 
higher outflow velocities are more likely to vary than those at lower velocities, and 
weaker BALs are more likely to vary than stronger BALs. The fractional change in 
BAL strength correlates inversely with the strength of the BAL feature, but does not 
correlate with the outflow velocity. Both the short-term and long-term data indicate 
the same trends. The observed behavior is most readily understood as a result of the 
movement of clouds across the continuum source. If the crossing speeds do not exceed 
the local Keplerian velocity, then the observed short-term variations imply that the 
absorbers are <6 pc from the central quasar. 

Key words: galaxies: active - quasars: general - quasars:absorption lines. 



1 INTRODUCTION 

Accretion disk outflows play a key role in the physics of 
quasars and their environs. Quasar outflows may be an in- 
tegral part of the accretion process and the growth of su- 
permassive black holes (SMBHs), by allowing the accret- 
ing material to release angular momentum. These outflows 
may also provide enough kinetic energy feedback to have 
an effect on star formation in the quasar host galaxies, 
to aid in 'unveiling' dust-enshrouded quasars, and to help 
distribute metal-rich gas to the interga l actic medium (e .g., 
iDi Matteo. Springel. fc HernguistlboOSl . IMoU et al.ll2007^ . 

The location and three-dimensional structure of quasar 
outflows are poorly understood. Sophisticated models have 
been developed that envision these outflows as arising from a 
rotating accretion disk, with acceleration to high speeds by 
radiative and/or magneto-centrifugal forces (jMurrav et al.l 
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1 19951 . IProea fc KallmanI l2004l . IProeal l2007l . lEverettl l2005l l. 
Improved observational constraints are needed to test these 
models and allow estimation of mass loss rates, kinetic en- 
ergy yields, and the role of quasar outflows in feedback to 
the surrounding environment. In this paper, we focus on 
the most prominent signatures of accretion disk outflows 
that appear in quasar spectra, the broad absorption lines 
(BALs). BALs are defined to have velocity widths >2000 
km s~^ at absorption depths >10% below the continuum 
ijWevmann et al.lll99ll), and they appear in t he spectra of 
~10-15% of quasars l|Reichard et al.l l2003al . [ Trump et al.l 
[2006) . Since only a fraction of quasar spectra display these 
features, the presence of BALs could represent a phase in the 
evolution of a quasar and/or particular orientations where 
the outflow lies between us and the quasar emission sources. 
Other studies have found support for the latter scenario, 
given the similarity of vari ous properties between BALQSOs 
and non- BALQSOs (e.g.. IWevmann et al.l[l99l] . IShen etlo] 
l2008h . 
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Studying the variability in these absorption features can 
provide important insight into the structure and dynamics 
of the outflows. For example, we can use information about 
short-term variability to place constraints on the distance of 
the absorbing material from the central SMBIf. The more 
quickly the absorption is varying, the closer the absorber is 
to the central source, b ased on nominally sh orter crossing 
times for moving clouds (|Hamann et al.ll2008l and §5 below) 
or the higher densities req uired for shorter recombination 
times (|Hamann et al.iri997l ). Long-term variability measure- 
ments provide information on the homogeneity and stabil- 
ity of the outflow. A lack of variability on long time-scales 
implies a smooth flow with a persistent structure. General 
variability results provide details on the size, kinematics, and 
internal makeup of sub-structures within the flows. Variabil- 
ity studies can also address the evolution of these outflows 
Eis the absorption lines are seen to come and go, or one type 
of ou tflow feature evolves into another (e.g.. iHamann et al.l 
I2OO8I ). 

A small number of previous studies have investigated 
variability in BALs, emphasizing Civ A1549 because of its 
prominence and ease of measurement. [Barlow ( 199S) car- 
ried out a study of 23 BALQSOs , covering time-scales of 
At <1 yfl iLundgren et al.l (|200it used Sloan Digital Sky 
Survey (SDSS) spectra of 29 q uasars to study Ciy BAL 
variability on similar time-scales. iGibson et al.l (|2008l ) stud- 
ied variability on longer time-scales (3-6 years) by combin- 
ing data for 13 BALQSOs from the Large Br ight Quasar 
Survey (LBQS) and SDSS. lGibson et al.1 (f201oh reports on 
variability on multi-month to multi-year time-scales, using 
3-4 epochs of data for 9 BALQSOs. They also compare vari- 
ability in Si IV absorption to variability in C iv absorption 
in their sample. 

In the present study, we go beyond existing work by 
carrying out a monitori ng program of BALQSOs from the 
sample of lBarlowl l|l993l ). This strategy provides BAL vari- 
ability data covering longer time-scales as well as a wider 
range in time-scales within individual sources. We include 
archival spectra from the SDSS, when available, to augment 
the temporal sampling. This paper reports our results for 
Civ a 1549 variability with measurements selected to en- 
able comparison between short-term (0.35 to 0.75 yr) and 
long-term (3.8 to 7.7 yr) behavior. We identify trends in 
the data with velocity and the depth of the absorption. We 
avoid using equivalent width (EW) measurements, which 
can minimize a change that occurs in a portion of a much 
wider trough. In subsequent papers, we will discuss variabil- 
ity properties in the entire data set, including more epochs 
on individual quasars, sampling in a much shorter time do- 
main, and comparisons between the Civ and Siiv A1400 
variabilities to place constraints on the outflow ionizations 
and column densities. Section 2 below discusses the obser- 
vations and the quasar sample. Section 3 describes the steps 
we followed to identify the C iv BALs and where they varied. 
Section 4 describes our results, and Section 5 discusses the 
implications of these results with comparisons to previous 
work. 



^ Throughout this paper, all time intervals are measured in years 
in the rest frame of the quasar. 



2 DATA AND QUASAR SAMPLE 

Betw een 1988 and 1992.lBarlowl l|l993t ) and his collaborators 
(e.g., iBarlow et al] 11993 7 obtained spectra for 28 BALQ- 
SOs with the Lick Observatory 3-m Shane Telescope. They 
obtained multi-epoch spectra for 23 of these quasars for a 
study of short-term BAL variability. They selected these ob- 
jects from already known BALQSOs, with redshifts of 1.2 
to 2.9. When selecting objects for their monitoring program, 
they gave preference to quasars known to have either absorp- 
tion line or photometric variability. However, in most cases, 
this information was not available, and they selected more 
quasars at higher redshifts due to higher detector sensitiv- 
ity at longer wavelengths. While this sample of BALQSOs 
is not randomly selected, it does cover a wide range of BAL 
stre ngths (s e e Tab le 1 and Figure 1). 

iBarlowl l|l993h used the Lick Observatory Kast spectro- 
graph to obtain spectra with settings for high resolution, 
R = A/AA ~ 1300 (230km s~^), and moderate resolution, 
R ^ 600 (530 km s~^). Most of the objects were observed 
at both settings. For our analysis of these data discussed 
below, we used the high-resolution observations for most 
sources. In some cases, only a moderate-resolution observa- 
tion is available or choosing the high-resolution observation 
would compromise wavelength coverage. Since BALs have a 
width of at least 2000 km s"'^, a resolution of 530 km s~^ is 
sufficient to detect changes in their profiles. 

Starting in January 2007, we have been using the MDM 
Observatory 2.4-m Hi ltner telescope to reobserve the BALQ- 
SOs in the sample of lBarlowl l| 19931 ). with the goal of moni- 
toring BAL variability over a wide range of time-scales, from 
<1 month to nearly 8 years. We also supplemented our data 
set with spectra from the Sloan Digital Sky Survey (SDSS). 
So far, we have over 120 spectra for 24 BALQSOs, and we 
continue to collect more data. We note that two of these 
objects are not strictly BALQSOs because they have a bal- 
nicity index of zero (see §3.2 below). Nonetheless, we include 
them in our sample because they do have broad absorption 
features at velocities that we include in our variability anal- 
ysis (§3.3). 

We used the MDM CCDS spectrograph with a 350 
groove per mm grating in first order and a 1" slit to yield 
a spectral resolution of f? « 1200 (250 km s~^), with wave- 
length coverage of ~ 1600 A. The spectrograph was rotated 
between exposures to maintain approximate alignment of 
the slit with the parallactic angle, to minimize wavelength- 
dependent slit losses. The wavelength range for each obser- 
vation was determined based on the redshift of the target 
quasar such that the coverage was blue enough to include 
the Lya emission and red enough to include the C iv emis- 
sion feature. The coadded exposure times were typically 1-2 
hours per source. 

We reduced the data using standard reduction tech- 
niques with the software. The data are flux calibrated 
on a relative scale to provide accurate spectral shapes. Ab- 



^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories 
(NOAO), which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under cooperative agree- 
ment with the National Science Foundation. 
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solute flux calibrations were generally not obtained due to 
weather or time constraints. 

The Sloan Digital Sky Survey is an imaging and spec- 
troscopic survey of the sky at optical wavelengths. The data 
were acquired by a dedicated 2.5-m telescope at Apache 
Point Observatory in New Mexico. The spectra cover the 
observer-frame optical and near- infrared, from 3800 to 9200 
A, and the resolution is R ~ 2000 (150 km s~^). For spectra, 
typically three exposures were taken for 15 minutes each, 
and more exposures were taken in poor conditions to achieve 
a target signal-t o- noise ratio. The multiple obser vations were 
then co-added l|Adelman-McCarthv et alll2008l ). The SDSS 
includes observations of 11 of the quasars in the Lick sample. 
However, 3 of the SDSS observations are unusable because 
the redshifts of those quasars are too low for the C iv BALs 
to appear in the SDSS spectra. Therefore, we have usable 
SDSS spectra for 8 of the objects in our sample, and they 
were taken between 2000 and 2006. 



3 ANALYSIS 

3.1 Short-Term and Long-Term Data Sets 

From the data in our BAL monitoring campaign, for each 
object, we select one pair of observations for the short-term 
analysis and one pair for the long-term analysis. For the 
long-term analysis, we want the longest time baseline possi- 
ble. The first observation for each object is the earliest Lick 
observation with the highest available resolution as well as 
wavelength coverage that extends from the Lya emission line 
through the C iv emission line. For the second observation 
for each object, we use the most recent MDM spectrum. We 
have long-term data for 23 objects, with resulting At of 3.8 
to 7.7 years. We note that for certain objects, the redshift is 
low enough that the Lya emission line is not shifted into the 
wavelength coverage of our data and therefore is not present 
in some of the spectra (e.g., 0946-1-3009). 

For the short-term analysis, the range in the values of 
At among the objects should be small, while including as 
many objects as possible. The best compromise we found 
was a range for At of 0.35 to 0.75 yr, which allows us to 
include 18 objects in the short-term analysis. In order to in- 
clude this many objects, we use data for some objects that 
only covers the region from the Si iv emission to the C iv 
emission, which is sufBcient for the analysis in this paper. 
These observations were selected separately from the long- 
term epochs and were chosen to fit into this At interval. 
These observations were taken before or at any time in be- 
tween the two long-term observations. There is one object 
for which we have short-term Lick data, but no long-term 
data (2225-0534), so we have in total 24 objects in our full 
sample. 



3.2 Characterizing the Quasar Sample 

Table 1 summarizes our full sample of long-term and short- 
term data, as described above. The first four columns pro- 
vide information on all 24 quasars. Zem is the emission red- 



shiflfl, and BI is the "balnicity index," as described below. 
The remaining columns are separated into the short-term 
analysis and long-term analysis. The columns with "vary?" 
have either Yes or No values to indicate variability anywhere 
in the C iv troughs (refer to §4), and At is the time difference 
(in years) in the quasar frame between observations. The ob- 
servations from 1988 to 1992 are from Lick, the observations 
from 2000 to 2006 are from SDSS, and the observations from 
2007 to 2008 are from MDM. 

Figure 1 contains spectra for all 24 objects, showing 
the long-term comparisons between a Lick spectrum and an 
MDM spectrum. The one exception is 2225-0534, for which 
we only have short-term Lick data. The velocity scale is 
based on the wavelength of C iv in the observed frame calcu- 
lated from the redshifts given in Table 1. We use a weighted 
average of the C iv doublet wave lengths, i.e., 1549.05 A . 

The balnicity index (BI, IWevmann et al.l Il99ll ) is 
designed to differentiate between BALQSOs and non- 
BALQSOs and provide a measure of the BAL strength, 
expressed as an equivalent width in units of velocity. It 
quantifies absorption that is blueshifted between —3000 and 
—25000 km s^^ from the Civ emission line centre, is at 
least ~10% below the continuum, and is wider than 2000 
km s~^. By this definition, if a quasar has a BI greater than 
zero, then it is classified as a BALQSO. The highest possible 
value of BI is 20,000. 

We measure BI in one observation per object to charac- 
terize the sample and allow for comparison with BALQSOs 
in other studies. We do not use BI in our variability analysis. 
In order to normalize the spectra and measure BI, and later 
absorption strength (see §3.3), we first had to fit a pseudo- 
continuum to one spectrum per quasar, which includes the 
quasar's continuum emission plus the broad emission lines. 
We adopt the Lick observations used in the long-term analy- 
sis as our fiducial epoch for which we fit the pseudo-continua 
and measure BI. We constructed this pseudo-continuum by 
fitting a power-law to the continuum away from any emis- 
sion or absorption lines, and then adding to that an esti- 
mate of the C IV and, if necessary, the Si IV broad emission 
hue (BEL) profilefl To fit the continuum, we fit a power- 
law function to two regions on the spectrum that are free of 
strong emission lines. The preferred wavelength ranges for 
the fits were 1270 to 1350 A and 1680 to 1800 A. These 
ranges were adjusted to avoid emission and absorption fea- 
tures as much as possible or due to the limits of the wave- 
length coverage. Fitting the C IV emission is complicated by, 
for example, potential asymmetries in the e mission profiles 
and ab sorption in the emission region (e.g.. iReichard et al.l 
l2003bl ). Therefore, we are fiexible in our fitting routine, us- 
ing between 1 and 3 gaussians to define the line profile. The 
gaussians have no physical meaning, and we simply use them 
for obtaining a reasonable fit to the profile. If there is any 
absorption on top of the emission, we ignore those wave- 
lengths when defining the fit. In cases where the absorption 
is wide enough that it obscures up to half of the emission 

The values of Zem were obtained from the NASA /IPAC Extra- 
galactic Database (NED), which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. 
* Note that the feature we attribute to Si IV may also include 
O iv] A1400 emission. 
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Figure 1. Smoothed spectra of all 24 quasars in our sample, showring the long-term comparisons between a Lick Observatory spectrum 
(bold curves) and the most recent MDM data (thin curves). The one exception is 2225-0534, for which we only have short-term data 
(see Table 1). The vertical flux scale applies to the Lick data, and the MDM spectrum has been scaled to match the Lick data in the 
continuum. The dashed curves show our pseudo-continuum flts. The thin horizontal bars indicate intervals of BAL absorption included 
in this study, and the thick horizontal bars indicate intervals of variation within the BALs. The formal la errors are shown near the 
bottom of each panel. 
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Figure 1. continued. 
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Figure 1. continued. 
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Table 1. Quasar Data and CIV Variability Results 
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These Lick observations were taken at the lower resolution setting (i? 600; see §2). 



line (e.g. 1413+1143), we assume the emission line is sym- 
metric. In cases where the gaussian, or multiple-gaussian, 
fit did not provide the right shape and we have enough in- 
formation about the profile (e.g. 1423+5000), we manually 
adjusted the fit to better match the C IV emission line. 

For quasars where the Civ absorption is at a velocity 
that overlaps the Si iv emission, we also fit the Si iv emission 
line. If most, or all, of the Siiv emission feature is absorbed, 
we synthesized a Silv profile by taking the Civ fit, shifting 
its centre wavelength to the location of Silv, adjusting the 
FWHM (given the wider separation between the Silv dou- 
blet lines), and scaling the emission strength by a factor of 
0.342. This sca le factor is based on comp osite quasar spec- 
tra from which FVanden Berk et all (|200ll ') found the typical 
strength of various emission lines. We relied on this scale 
factor in cases such as 0146+0142, where we have little in- 
formation about the Silv emission line. Otherwise, in a case 
such as 0846+1540, we fit between 1 and 3 gaussians to the 
Si IV line, as we did for the Civ emission, because we have 
more information about the profile of the Si iv emission line. 
The pseudo-continuum fits for each object are shown by the 
smooth dotted curves in Figure 1. 



3.3 Measuring the BAL Variability 

For each quasar, we compare the two spectra in our short- 
term and long-term intervals (see Table 1) to search for vari- 
ability in the C iv BALs. First, we determined the veloc- 
ity ranges over which BAL absorption occurs. These ranges 
might differ between the two observing epochs. We therefore 



set the range to the extremes observed in either epoch. Our 
definition of these ranges is guided by the definition of BI, 
i.e. they must present contiguous absorption reaching ^10% 
below the continuum across >2000 km s-\ We do not follow 
the requirement that the absorption be blueshifted between 
—3000 and —25000 km s~^ because we want to include all 
Civ broad absorption in our analysis (as a result, our sam- 
ple contains two objects, 0302+1705 and 0846+1540, with a 
BI of 0). In general, the Civ BALs occur between the Silv 
and C IV emission lines, but some of the C iv absorption is 
blue-shifted enough that it appears blueward of the Silv 
emission line center (e.g., 0146+0142). In these cases, we 
can confirm that the absorption blueward of the Silv emis- 
sion line is due to C iv, and not Si iv, because any significant 
Si IV absorption should be accompanied by C iv absorption 
at the same velocity. This is evident from numerous observa- 
tions of BALQSOs where the Silv and Civ absorption are 
clearly identified (e.g.. Figure 1, Korista et al. 1993), and 
from theoretical considerations of the relative strengths of 
these lines (Hamann et al. 2008). Our sample includes broad 
absorption troughs from to —40000 km s~^. 

In order to compare any two spectra for variability in 
the defined absorption ranges, we scaled the spectra so that 
regions free of absorption and emission in each spectrum 
match. This is necessary because of uncertainties in the flux 
calibration and possible real changes in the quasar emissions. 
We applied a simple vertical scaling so that the spectra from 
each epoch match the fiducial Lick spectrum along the con- 
tinuum redward of the Civ emission feature (i.e., from 1560 
A to the limit of the wavelength coverage), between the Si iv 
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and C IV emission ('^1425-1515 A), and between the Lya and 
the Siiv emission (^1305-1315 A). In nearly all cases, a sim- 
ple scaling produced a good match between the spectra from 
different epochs. The results for the two long-term epochs 
are shown in Figure 1. In some cases, however, there were 
disparities in the overall spectral shape of the two spectra. 
To remove these disparities, we fit either a linear function 
(for 0903-f 1734, 1413-f 1143, and 1423+5000) or quadratic 
function (for 00194-0107 and 1309-0536) to the spectral re- 
gions that avoid the BALs in a ratio of the two spectra. We 
then multiplied this function by the comparison spectrum. 

We identify BAL variability in velocity intervals that 
correspond to at least 2 resolution elements in the lower res- 
olution Lick spectra, i.e., 1200 km s~^. Our procedure to 
identify variability intervals relies mainly on visual inspec- 
tion because photon statistics alone are not sufficient for 
defining real variability. Most of the occurrences of variabil- 
ity recorded for this study are obvious by simple inspection. 
However, we also examined all cases of weak or possible vari- 
ability to ensure that no occurrences of significant real vari- 
ability are missed. We determined the significance of vari- 
ability based on two criteria. First, we calculate the average 
ffux and associated error within each variability interval and 
use that to place an error on the ffux differences between the 
two spectra. We include in this study only the intervals of 
variability where the fiux differences are at least 4-sigma. 
We emphasize that all of the intervals which varied by 4- 
sigma or more are readily identified by our initial inspection 
technique. Therefore, there is no danger of variable intervals 
being missed. However, the 4-sigma threshold alone can in- 
clude spurious or highly uncertain variability results because 
of uncertainties not captured by the photon statistics. These 
additional uncertainties can arise from the flux calibration, a 
poorly constrained continuum placement, underlying emis- 
sion line variability or our reconstruction of a severely ab- 
sorbed emission line profile. We therefore apply a second 
significance criterion based on our assessment of these ad- 
ditional uncertainties. We take a conservative approach by 
excluding variability intervals where the additional uncer- 
tainties might be large (even if the flux difference passes the 
formal 4-sigma threshold described above). This second cri- 
terion caused the exclusion of just a few intervals. See below 
for examples. The velocity intervals finally classified as vari- 
able in the long-term data are indicated by bold horizontal 
lines in Figure 1. 

To help clarify our assessments of the additional uncer- 
tainties, we point out some specific velocity intervals that 
pass the 4-sigma threshold and appear variable by casual 
inspection, but are nonetheless excluded from our analysis 
because the additional uncertainties are significant. For ex- 
ample, in 1246-0542, we do not include the region around 
— 17,000 km s~^ because the noise level in that portion of 
the trough makes it a very tenuous candidate for inclusion in 
the study. The region around -21,000 km s"^ in 1435-1-5005 
is a similar case because here the increased noise level of 
the spectra towards greater outflow velocities makes a de- 
termination of variability more uncertain. We note that the 
increased scatter in the data for these two objects is less 
noticeable in the smoothed spectra plotted in Figure 1. In 
a couple cases, a slight difference in the slope between two 
comparison spectra adds some uncertainty. In 1413+1143, 
the blue wing of the C iv trough in the MDM spectrum ex- 
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Figure 2. The average normalized absorption strength {(A)), cal- 
culated throughout all of the BALs in the quasars in our long-term 
subsample, as a function of the outflow velocity of the absorbing 
material. The individual points in each bin (with width 1000 km 
s~^) represent the {A) for a single quasar, and the open squares 
represent the average of all points in that bin. 



tends slightly above the deflned continuum. This may be due 
to some anomaly in the fiux calibration, so we define the blue 
end of the BAL (and the variability interval) as —20,600 km 
s~^. In 2225-0534, there is very little continuum available for 
matching the two comparison spectra, and this added uncer- 
tainty makes the small fiux difference at the bluest portion 
of the BAL trough too tenuous to be included. A slight dif- 
ference in slope between the two epochs could be responsible 
for any apparent changes in the wing of this BAL, and we 
do not have enough visible continuum to discount this pos- 
sibility. 

Identifying variability in the spectral regions which 
overlap the C iv and Si iv emission lines is further com- 
plicated by possible variability in the emission lines them- 
selves fe.g.. lWilhite et a"Lll2005l ). We record BAL variability 
in these spectral regions only if it is large compared to any 
reasonable change in the broad emission lines and/or the ab- 
sorption changes abruptly over just portions of a BEL pro- 
file. For example, in 0146+0142, the interval from -38,800 to 
—32,700 km s^^ clearly shows BAL absorption weakening to 
approximately continuum level. In another case, 0932+5006, 
we include the intervals from —37,500 to —35,400 km s~^ 
and —28,400 to —27,000 km s"'^ because there is clearly 
BAL absorption there in the MDM spectrum. We did not 
attribute the change in fiux from approximately —34,000 to 
-30,000 km s"^ to a BAL because it is unclear what is the 
primary cause of variability in this interval. While the iden- 
tifications of BAL absorption and variability in the regions 
overlapping C iv and Si iv emission are secure, the measure- 
ments of absorption strength {{A)\ see below) in this re- 
gion have added uncertainty since fitting the emission lines 
is more uncertain than fitting the continuum. When look- 
ing at correlations with velocity in this study, we identify 
trends using all the data, but also confirm the existence of 
these trends in the velocity range —27,000 to —8,000 km 
s~^, which avoids the emission lines. For correlations with 
absorption strength (see below), since there are few cases of 
variability at the extremes of the examined velocity range, 
their inclusion does not significantly alter the plots or the 
final results. 

In order to look for trends in the data, we calculated the 
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absorption strength as a function of velocity in each spec- 
trum. We first normahzed each spectrum using the pseudo- 
continuum fit (see §3.2) for that quasar. We can apply the 
pseudo-continua derived in §3.2 to the other epochs for each 
quasar because those epochs have already been scaled to 
match the fiducial Lick spectrum. However, we adjust the 
emission line fits in cases where the emission line varied. We 
define absorption strength, A, as the fraction of the nor- 
malized continuum flux removed by absorption (0 ^ A ^ 
1) within a specified velocity interval. A single A value is 
adopted for each velocity interval based on the average flux 
within that interval. If there is variability, then the absorp- 
tion strength changed and we calculate the average absorp- 
tion strength, (A), and the change in absorption strength, 
AA, between the two epochs being compared. The veloc- 
ity intervals used to calculate {A) and AA are deflned to 
i) include only portions of the spectrum with BAL absorp- 
tion and ii) clearly separate the spectral regions within BAL 
troughs that did and did not vary. Specifically, the variable 
and non-variable absorption regions are each divided into 
equal-sized bins, with a bin width of 1000 to 2000 km s~^ 
depending on the velocity width of the region. This strat- 
egy makes the bin widths used to measure (A) and AA 
similar to the nominal width of 1200 km s~^ used above, 
while providing the fiexibility needed to begin and end at 
the boundaries defined above between absorbed/unabsorbed 
and variable/non- variable spectral regions. In each of these 
velocity bins, we also calculate the fractional change in the 
absorption strength, |Aj4|/(A), which can have values from 
to 2. Spectral regions that did not vary according to the 
criteria described above have AA set to zero automatically. 
Note that, with (A) and AA determined in this way, each 
quasar can contribute more than once to these quantities in 
our variability analysis, depending on the range of strengths 
and velocities covered by its BAL troughs, because each bin 
showing variability is treated as a separate occurrence. 

We plot the relationship between absorption strength, 
(A), and outfiow velocity for all the quasars in our long- 
term subsample in Figure 2. A correlation between these 
two quantities will have implications for our results when 
we compare variability measures to each of these quantities 
individually (see §4 and §5). There is a trend for lower val- 
ues of (A) at increasing outflow velocity, which is especially 
apparent from —27,000 to —8,000 km s~^, away from the re- 
gions of Si IV and C iv emission. The errors in the individual 
points are much smaller (cr^^^ ^0.01) than the dispersion of 
points in each bin. However, the bins at <-30,000 km s^^ 
only contain one or two quasars, so the average (A) in those 
bins can be uncertain by up to a factor of 2, given the ad- 
ditional uncertainty of fltting the Siiv emission line. The 
weaker absorption found at higher velocities in our sample 
is consisten t with absorpt i on in the mean BAL spectra con- 
structed bv lKorista et al.l l|l993h . 



4 RESULTS 

Table 1 indicates whether there was significant variability in 
the C IV BALs of each quasar in the short-term and long- 
term data. We found that 39% (7/18) of quasars with short- 
term data varied, whereas in the long-term, 65% (15/23) 
varied. One quasar varied only in the short-term, so com- 
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Figure 3. The top two panels show the number of quasars with 
C IV BAL absorption and the number with C IV BAL variability at 
each velocity. The third panel is the second panel divided by the 
top one, which gives the fraction of C IV BALs that varied at each 
velocity. The dashed, dotted, and solid curves indicate results for 
the long-term, short-term, and combined data set comparisons, 
respectively. The bottom panel displays the combined data set 
with la error bars. 



bining both short-term and long-term data, 67% (16/24) of 
the quasars varied overall. 

We investigated C iv BAL variability as a function of 
outflow velocity. Using the velocity ranges for each quasar 
where we have identified absorption and variability (§3.3), 
we count how many quasars have BAL absorption and vari- 
ability in each velocity bin. In Figure 3, the top two panels 
show the number of quasars with C iv BAL absorption and 
the number with C iv BAL variability at each velocity. The 
third panel is the second panel divided by the top one, which 
gives the fraction of C iv BALs that varied at each velocity. 
The dashed, dotted and solid curves indicate results for the 
long-term, short-term and combined data set comparisons, 
respectively. The plot implies that BALs at higher velocities 
are more likely to vary than those at lower velocities. 

The bottom panel in F igure 3 displays la err or bars 
based on lWilsonI (|l927l ) and lAgresti fc Coulll l|l998h . These 
errors are based on counting statistics for the number of 
quasars with absorption and variability at each velocity. We 
perform a least-squares fit to confirm that a correlation ex- 
ists between incidence of variability and velocity. The slope 
of the combined data is —1.51 ± 0.25 x 10"'', in the formal 
unit, fraction per km s~^. The slope is non-zero at a 6-sigma 
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significance. If we repeat tlie least-squares fit witli only the 
bins witli tlie most data (bins containing at least 8 quasars 
with absorption), the slope is — 1.76±0.33x 10~^ fraction per 
km s~^, which is non-zero at a 5-sigma significance. Using 
this restriction limits the fit to only the data from —29,000 
to —2,000 km s~^. We perform the least-squares fit a third 
time, where we restrict the fit to only the data in the range 
—27,000 to —8,000 km s^^ that completely avoids overlap 
with the Civ and Siiv broad emission lines. In this case, 
the resulting slope is —1.46 ± 0.54 x 10^'' fraction per km 
s~^, which is non-zero at a 2.7-sigma significance. Therefore, 
this trend is evident even when we remove the data at the 
lowest and highest velocities, where additional uncertainty 
from the BELs can affect the results. 

An important question is whether the trend seen in Fig- 
ure 3 is a correlation with velocity in an absolute sense, or 
the cumulative result of differences between the red versus 
blue-side behaviors in individual absorption troughs. To ad- 
dress this question, we divided the velocity range of absorp- 
tion for each quasar into a red and a blue half and plotted 
the incidence of variability separately for the red and blue 
portions. The result in both cases is consistent with the be- 
havior shown in Figure 3. We do not sec a higher incidence 
of variability in the blue portions of troughs versus the red 
portions. This indicates that the trend in Figure 3 is a trend 
in the ensemble dependent on the absolute velocity rather 
than the relative velocity within a given trough. Consistent 
with this finding is the fact that none of the BALs varied at 
a velocity between —3500 km s~^ and km s~^, while the 
two BALs in our sample at <— 33,000 km s~^ both varied. 

Another factor that might affect the results in Figure 
3 is that BALs at higher velocities tend to be weaker (Fig- 
ure 2). In Figure 4, we examine Civ BAL variability as a 
function of absorption strength. Using the values of (A) cal- 
culated in §3, and a bin size of 0.05, we count how many 
quasars have absorption in each bin. We then count how 
many quasars have at least one occurrence of variability in 
each absorption strength bin. Each quasar is counted at most 
once in each bin, even if it has the same absorption strength 
at more than one location in the spectrum. As in Figure 3, 
the third panel is the second panel divided by the top one, 
in this case yielding the fraction of quasars with a BAL that 
varied in each absorption strength bin. Figure 4 indicates 
that quasars are more likely to exhibit BAL variability at 
weaker absorption strengths. As in Figure 3, we calculate 
errors for each absorption strength value. We again perform 
a least-squares fit, which gives a slope of —0.487 ± 0.095 
fraction per unit absorption strength. The slope is non-zero 
at a 5-sigma significance. We repeat the least-squares fit 
with only the bins with the most data, as for Figure 3; this 
restricts the fit to all of the bins except the one at {A) > 
0.95. The resulting slope is —0.476 ± 0.100 fraction per unit 
absorption strength, which is non-zero at a 4.8-sigma signif- 
icance. 

Figure 4, however, does not account for multiple occur- 
rences of the same absorption strength within an individual 
quasar. In §3.3, we define a method of dividing each BAL 
into ~1200 km bins, and we now consider these bins 
as individual occurrences of absorption. In Figure 5, wc plot 
the number of these occurrences at each absorption strength 
value in the top panel and then the number of these occur- 
rences that varied in the second panel. The third panel is the 
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Figure 4. Same as Figure 3 except the quantities are plotted as a 
function of the normalized absorption strength, averaged between 
the epochs being compared. 

second panel divided by the top one, showing the fraction 
of occurrences at each absorption strength that varied. An 
individual BAL may be counted more than once at each ab- 
sorption strength value because we consider portions of each 
BAL trough in this figure. The bottom panel is the same 
as the third panel, except it shows just the long-term data 
with error bars, calculated as in Figures 3 and 4. A least- 
squares fit for the long-term data in Figure 5 gives a slope of 
—0.380 ± 0.052 fraction per unit absorption strength, which 
is non-zero at a 7-sigma significance. Figure 5 confirms that 
weaker portions of BAL troughs are more likely to vary than 
stronger ones. 

Figures 3 to 5 compare the incidence of variability to 
properties of the BALs. Next, we use the values of | Aj4|/(yl), 
calculated in §3, to compare the amplitude of variability to 
the same BAL properties. 

In Figure 6, we plot this fractional change in absorp- 
tion strength as a function of the outflow velocity, using bin 
sizes of 1000 km s~^ width, as in Figure 3. In each velocity 
interval, the individual points represent variability in a dif- 
ferent quasar. The average values of lAyl|/{A) for just the 
variable (non-zero) points in each bin are shown by the open 
squares. The asterisks show the average if the non-varying 
absorption intervals are included. This plot does not indicate 
a statistically significant trend between fractional change in 
absorption strength and the outfiow velocity. 

In Figure 7, we plot the fractional change as a function 
of average normalized absorption strength, similar to Figure 
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Figure 5. The top two panels show the total number of occur- 
rences of absorption, and the number of those occurrences that 
varied, at each strength. The third panel is the second panel di- 
vided by the top one. The bottom panel displays the long-term 
data with Icr error bars. 
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Figure 6. Fractional change in the depth of BAL absorption as 
a function of velocity. The symbols follow the same pattern as 
Figure 2, with individual points representing variability in a sin- 
gle quasar. Here, the open squares indicate the average of these 
variable, non-zero points in each velocity bin. The asterisks show 
the average if the non- varying absorption intervals are included. 
The top panel shows the short-term results, while the bottom 
panel plots long-term data. We note that the maximum value of 
Ayl|/(A) is 2 (see Figure 7), which indicates absorption appear- 
ing or disappearing completely. 



6. In this plot, however, individual quasars will contribute 
multiple times to the measurements at a given absorption 
strength if the quasar has variability in both the red and blue 
sides of an absorption trough and/or in multiple troughs. 
The individual points represent variability at different ve- 
locities in each quasar. The average of all the variable (non- 
zero) points in each bin are shown by the open squares, while 
the asterisks show the average |Aj4|/{A) in each bin if non- 
varying absorption regions are included. This plot indicates 
that the fractional change in strength is greater in weaker 
features. However, the errors in the values of |Ay4|/{yl) are 
correlated with the values of {A). At an {A) of 0.15, the 
typical error values range from ~0. 11-0.33, whereas at {A) 
of 0.60, the error values are ~0. 02-0. 06. A larger error at 
smaller values of {A) should cause a larger dispersion in 
the points in that region of Figure 7, but instead the val- 
ues of |Ay4|/{j4) converge towards the maximum possible 
value of |AA|/(^). The dispersion of points due to errors in 
|A^|/(^) alone cannot explain the trend in this figure. In- 
stead, this figure indicates that we see weaker absorbers, or 
weaker portions of BALs, appear or disappear completely, 
but strong features do not typically appear or disappear over 
the time-scales that our data cover. 

Beyond looking at trends with velocity and absorption 
strength, it is instructive to examine the relationship be- 
tween these two parameters. Figure 2 displays {A) as a 



function of velocity, revealing that these quantities are cor- 
related. The lack of a clear trend in Figure 6 is surprising 
given that weaker systems are more variable (Figure 7) and 
they tend to appear more often at higher velocities (Figure 
2). Evidently, the correlation of |AA|/{j4) with absorption 
strength is stronger than the correlation with velocity. We 
also investigate how |Ayl| varies with velocity and with {A). 
We do not include the plots here, but we find no trend be- 
tween |Ayl| and velocity or absorption strength. 

We do not find any significant differences in the trends 
described above between the short-term and long-term data. 
The main difference is in the incidence of variability, 39% in 
the short-term versus 65% in the long-term data. The long- 
term data also show slightly larger changes in absorption 
strength, i.e., |AA|, compared to the short-term data. The 
average value of |Ayl| is 0.16 ± 0.08 in the short-term and 
0.22 ± 0.10 in the long-term. 

We also note again that we include two quasars 
(0302-f 1705 and 0846-1-1540) from the Lick sample l|Barlowl 
119931 ) that are not technically BALs. However, these two 
objects do not significantly affect the results. If we omitted 
them, then the fraction of quasars with variability would 
be 35% (6/17) in the short-term and 67% (14/21) in the 
long-term. In the case of 0302-1-1705, it does not qualify as a 
BALQSO because its broad absorption feature is at low ve- 
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Figure 7. Same as Figure 6, except tlie fractional cliange is now 
plotted as a function of average normalized absorption strength. 
In this plot, individual quasars can contribute multiple times to 
the measurements at a given absorption strength. The solid curve 
represents the maximum possible value of |AA|/(j4). 

locities (centered at r--^— 2100 km s~^). It also has {A) reach- 
ing >0.90. It is a strong feature at low velocity that does not 
vary, fitting the trends discussed above. In 0846+1540, there 
is broad absorption at low velocity (centered at ~— 2900 km 
s~^), with average {A) of r--^0.27, and at high velocity (cen- 
tered at ~-26,800 km s~^), with average {A) of ~0.18. This 
also fits the overall trend for higher incidence of variability 
at higher velocities. 



5 SUMMARY AND DISCUSSION 

We have analyzed short-term and long-term C iv A 1549 BAL 
variability in a sample of 24 quasars, and we looked for 
trends with the outflow velocity and absorption strength. 
We found that 39% (7/18) of the quasars varied in the 
short-term, intervals of 0.35 to 0.75 yr, whereas 65% (15/23) 
varied in the long-term, intervals of 3.8 to 7.7 yr. Overall, 
67% (16/24) of the quasars varied. We find a trend for vari- 
ability to occur more often at higher velocities (Figure 3) 
and in shallower absorption troughs (or shallower portions 
of absorption troughs; Figures 4 and 5). When looking at 
the fractional change in strength of the varying absorption 
features, there is no apparent significant correlation with ve- 
locity (Figure 6), but there is a trend for a larger fractional 
change in absorption strength in shallower features (Figure 
7). We do include two objects with BI=0, but the variability 
in the broad absorption in these two quasars is consistent 
with the trends described above. 



iGibson et"al] (|2008l ') studied Civ BAL variability on 
time-scales of At~3-6 yrs in 13 quasars with two epochs 
of data and found that 92% (12/13) varied. This is larger 
than our overall percentage of 67%. Th is is sur p rising be- 
cause the selection criteria adopted by iBarlowl (|l993l ) for 
the sample used here might be biased toward more variable 
sources (§2). Evidently, this bias is small or nonexiste nt. The 
difference between our result and Gibs on et al. 1 12008') might 
be due to the small numbers of objects studied and/or the 
more conservative approach we took for identifying variable 
abs orption. 

iLundgren et al.l l|2007t ). which is a study of BAL vari- 
ability on time-scales of <1 yr, find an inverse correlation 
of fractional change in equivalent width (EW) with average 
EW. While EW is a different calculation than our measure 
of absorption strength, {A), encapsulating the entire BAL 
profile in one number instead of considering separately the 
different profile regions (§3.3), this result is consistent with 
our findings that fra c tiona l change in {A) correlates with 
{A). iLundgren erakl l|2007l ) also find no overall trend be- 
tween fractional change in EW and outflow velocity, which 
is again consistent with our results. While not plotted here, 
we do not find a significant correlation between |Ay4| and 
velocity or absorption str e ngth, w hich is consist e nt wit h the 
analysis of ICibson et all l|2008l 'l. ICibson et akl (|2008l 'l also 
calculate the change in absorption strength, with values con- 
centrated in the range 0.15-0.25. We find an average value 
of j A^l fo r the long-term data of 0. 22 ± 0.10. 

Both I Gibson et al] (|2008h and ILundgren erall (|2007l ) 
comment on the lack of evidence for acceleration in the BAL 
features they study. We also find no evidence for accelera- 
tion in our sample. The BAL variations only involve the ab- 
sorption trough growing deeper or becoming shallower. We 
note, however, that limits on the acceleration are difficult 
to quantify for BALs because, unlike narrower absorption 
lines, they do not usually have sharp features to provide an 
accurate velocity reference. BALs also exhibit profile vari- 
ability, which can shift the line centroid without necessarily 
corresponding to a real shift in the velocity of material in 
the outflow. 

In this work, we compare two different time intervals, 
and we do not find a significant difference between the short- 
term and long-term data in the trends described in §4. How- 
ever, the incidence of variability and the typical change in 
strength is greater in the long-term than in the short-term. 
This indicates that BALs generally experience a gradual 
change in strength over multi-year time scales, as opposed 
to many rapid changes on time-scales of less than a year. 
This will be investigated more in a later paper where we 
will include all the epochs from our data set. 

Variability has been investigated in other types 

of absorbers th at are po tentially r e lated to BALs. 

[Narayanan et al.l (|2004l ') and IWise et all (|2004h looked for 
variability (At~0.3-6 yrs) in low-velocity narrow absorp- 
tion lines (NALs) and found that 25% (2/8) and 21% 
(4/19), respectively, of t he N ALs they studied varied. 
iRodngucz Hid algo et al] |2010h examined variability in 
mini-BALs, which are absorption lines with widths larger 
than those of NALs and smaller than those of BALs. They 
observed Civ mini-BALs in 26 quasars on time-scales of ~1 
to 3 years and found variability in 57% of the sources. This 
number is similar to our result for BALs, however this com- 
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parison is complicated by a differing number of epochs and 
diff erent time baselines. 

iRodrfguez Hidalgo et al] (|2010l ) found no correlation 
between incidence of mini-BAL variability and absorption 
strength or outflow velocity. However, this comparison might 
be skewed by the fact that mini-BALs tend to be weaker 
th an BALs. For example, none of the mini-BALs studied 
by iRodn'guez Hidalgo et all l|2010 !') have depths A > 0.6, 
while some of the BAL features in our sample have yl ~ 1. 
There are also differences i n the velocity distributio n s. The 
mini-BAL distribution in IRodrfguez Hidalgo et all (|201(]| ') 
peaks at ~— 20,000 km s~^, whereas the distribution of BAL 
absorption in our data peaks at —15,000 to —12,000 km 
s~^ (Figure 3). More work is needed, e.g., with larger sam- 
ples, to control for these differences in the absorption char- 
acteristics and compare BAL and mini-BAL variabilities on 
equal footing. Comparisons like that could provide valuable 
constraints on the physical similarities and relationship be- 
tween BAL and mini-BAL outflows. 

Two possibilities for the cause of BAL variability are 
the movement of gas across our line of sight and changes 
in ionization. We found a trend for lines at higher veloc- 
ity to vary more often, which suggests movements of clouds 
since faster moving gas might vary more. Moreover, this is 
a global trend, dependent on velocity in an absolute sense, 
rather than relative velocity within a given trough. However, 
it is possible that the root cause of the variability is associ- 
ated with the strength of the lines, and weake r lines happen 
to ap pear at higher velocities (see Figure 2 and lKorista et ahl 
Il993f ). In either case, faster moving material tends to have 
lower optical depths or cover less of the continuum source. 
Measured values of {A) thus provide either a direct measure 
of optical depth, or covering fraction in the case of satu- 
rated absorption. With the current results, we cannot yet 
disentangle the trends in velocity and {A) with certainty. 

Evidence from BAL variability studies nonetheless sup- 
ports an interpretation of varying covering factor. In ad- 
dition to greater variability at higher velocities, we found 
changes over small portions of the BAL troughs, rather than 
entire BALs varying, which can be understood in terms of 
movements of sub-structures in the flows. Changes in the 
ionizing flux should globally change the ionization in the 
flow, and we would expect these changes to be apparent over 
large portions of the BAL proflles an d not in just some small 
velocity range. iHamann et al.l (|2008l ) looked at the variation 
in the C iv and Si iv BALs in one quasar and found that the 
BALs were locked at roughly a constant Siiv/Civ strength 
ratio. They attribute this result to the movements of clouds 
that are optically thick in both lines. If the level of ioniza- 
tion is constant between observations or if r S> 1, then the 
changes in absorption strength in our sample are directly 
indicative of changes in the covering fraction. We will exam- 
ine the Silv/Clv line ratio in our BAL data in a subsequent 
paper. 

On the other hand, studies of NALs support the pos- 
sibi lity of changes in i oniz ation causing the line variabil- 
ity. [Mis^ra^eral] (120071 ) and lHamann et al.1 (|2010l ') observed 
qucisars where multiple NALs varied in concert, which sug- 
gests global changes in the ionization of the outflowing gas. 
If a change in ionization is causing the variability in the 
BAL outflows, then the optical depth (r) of the BALs is 
changing. Features with lower optical depth are much more 



sensitive to changes in ionization than those with higher 
optical depth. The outflowing gas is presumably photoion- 
ized by the flux from the continuum source, and changes 
in the continuum flux could cause a change in the ioniza- 
tion of the outflowing gas. However, studies to date have 
not found a strong corr elation between continuum variabil- 
ity and BAL variab il itv jLundgren et al.ll2007l . ICibson et al.l 
"iool, iBarlowl 1 1991 iBarlow et al.n i992^. casting doubt on 
ionization changes causing BAL variability. 

Finally, we use the bolometric luminosities of the 
quasars in our sample to derive characteristic time-scales 
of the flows for comparison to the observed time-scales. The 
bolometric luminosities for the quasars in our sample range 
from ~ 2 X 10*® to 3 X 10*^ ergs s~^, and the average is 
~ 7 X W^^ ergs s"'^. These values are based on the observed 
flux at 1 450 A ( rest-frame) in absolute flux-calibrated spec- 
tra from iBarlow (1993), a cosmology with Ho ~ 71 km s^^ 
Mpc, Qm ~ 0.3, Ha = 0.7, and a standard bolometric cor- 
rection factor L « 4.4ALa(1450A). Based on the average 
value of L ~ 7 X W^^ ergs s~^, a characteristic diameter 
for the continuum region at 1550 A is -D1550 ~ 0.006 pc 
and for the C iv BEL region is Dciv ^ 0.3 pc, assuming 
L = l/3Ledd and Mbh = 1.4x10^ M^ta llPeterson et al.ll20o4 . 
iBentz et al.ll2007l . iGaskelllliooil . iHamann fc Simonll201Gl ). If 
the outflowing gas is located just beyond the C iv BEL re- 
gion and it has a transverse speed equal to the Keplerian 
disk rotation speed, then the gas cloud would cross the en- 
tire continuum source in ~1 yr. The typical time-scale in our 
short-term data is ~0.5 yr and the typical change in absorp- 
tion strength is at least 10%, which nominally requires the 
moving clouds to cross at least 10% of the continuum source. 
This implies transverse speeds >1000 km s~^ and thus ra- 
dial distances that are conservatively <6 pc if the transverse 
speeds do not exceed the Keplerian speed. A typical time- 
scale in the long-term data is ~6 yr, which is enough time to 
cross the continuum source ~7 times (if the gas is located 
just beyond the Civ BEL region). This suggests a fairly 
homogeneous flow for at least the cases where there is no 
long-term variability. 

To place better constraints on the outflow properties 
and understand the underlying cause(s) of the line variabil- 
ities, our following papers will make more complete use of 
our entire data sample. We will examine the ratio of Siiv 
to C IV absorption to gain better insight into the cause(s) 
of the variability, and we will utilize more complete time- 
sampling instead of just looking at two specific ranges of 
time intervals. We have also obtained more data on BAL 
variability at the shortest time-scales. At ~ 1 week to 1 
month, in order to determine whether there is a minimum 
time-scale over which variability occurs. This will provide a 
better constraint on the location of the outflowing gas. 
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